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PREFACE 


The present thesis an out come of researches carried out by me 
in the field of “STUDY OF RADIATIVE & ELECTROMAGNETIC 
DISCONTINUTIES” under the supervision of Dr. S. P. Singh Bhadauria, M.Sc., 
Ph.D. Reader & Head, Deptt. of Maths, Pt. J. N. College, Banda (U.P.) affiliated 
to Bundelkhand University, Jhansi, is being submitted for the award of Ph.D. 
degree in Mathematics. This thesis has been devided into seven chapters, 
each chapter has been subdevided into a number of articles. 

The first chapter is introductory. It gives in brief, an idea about 
shock waves, Lagrangian Method, Eulerian Method, equation of state, equation 
of continuity, shock wave their existance, spherical and cylindrical shock wave, 
equation of motion and jump condition in magnetohydrodynamics, radiation 
phenomenon, radiant pressure, radiation flux, similarity principle and self 
similar gas motion and concept of self gravitation. 

The second chapter is devoted to unify study of PROPAGATION : 
OF PLANE RADIATIVE SHOCK WAVE IN UNSTABLE HOMOGENOUS ; 
MEDIUM taking disturbance headed by an radioactive shock surface of a 
rotating thunder in the earth medium and definded study is made for plane 

shock wave. | 

i 

t 

The Third chapter deals with ANALYTICAL SOLUTION OF SELF 
GRAVITATING GAS BEHIND A CYLINDRICAL SHOCK WAVE IN 
MAGNETOGASDYNAMICS. The motion of gas is assumed under azimuthal 
magnetic field, with gas in self gravitating. Energy of wave increases with time. 
An analytical study is made for cylindrical shock wave with azimuthal magnetic 
field. 



!n fourth chapter, an attempt is made for the ANALYTICAL 
SOLUTION OF SPHERICAL SHOCK WAVE IN A ROTATING GAS WITH 
AXIAL COMPONENT OF MAGNETIC FIELD, taking density as constant ahead 
the shock front, By using C.C.W. method, analytical solution are have been 
obtained for shock velocity and shock strength for weak and strong magnetic 
field. For strong shock, also we have considered two cases i.e. when the 
magnetic field is strong and when weak i.e. nonmagnetic case. 

Fifth chapter deals with ANALYSIS OF MAGNETO 
HYDRODYNAMICS SPHERICALLY SHOCK WAVE IN A MEDIUM WITH 
VARIABLE DENSITY, energy is suddenly released with spherical symmetry in 
a conducting gas, under the effect of azimuthal magnetic field where density 
varies with time. Numerical integration of differential equation is carried out in 
DEC 1090, Computer System at IIT, Kanpur by well known R.K.G’s program 

In sixth chapter, an attempt is made for the study of 
PROPAGATION OF PLANE SHOCK WAVE IN MAGNETO GAS DYNAMICS 
under isothermal condition where temperature gradient is zero and radiation 
effect have been taken out, Numerical integration of differential equation is 
carried out on DEC 1090, Computer at IIT Kanpur, by R.K.G’s program. 

In seventh chapter, an attempt is made for the analytical solution 
of RADIATIVE EFFECT OF POINT EXPLOSION IN STELLAR BODIES. Since 
the temperature of material is very high the radiation effect can not be ignore. 
The wave is produce by sudden point explosion. The effect of transverse 
magnetic field be taken out. 
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CHAPTER - 1 


INTRODUCTION 

The study of 'Fluid Dynamics" is directed to the behavior of a fluid 
in motion. The liquid and gas states are referred to generally as fluids. Some of 
tliG nolnblc oxamplos of application ate : flow of water through pipes; motion of 
an aircrafts or a missiles in the atmosphere. The study also yields the methods 
and devices for the measurement of various parameters, e.g. the pressure and 
velocity in a fluid at rest or in motion. 

The first notable works on fluids appeared in the seventeenth 
century. In 1687, Newton in his book ‘PRINCIPIA’ deals with the Influx of a 
fluid, the resistance of a fluid and the resistance of projected bodies and the 
wave motion. Daniel Bernoulli in 1738 determined the relation between 
pressure and velocity and formalized this in his theorem. In 1743 his son John 
Bernoulli applied the momentum principle to infinitesimal elements. Euler 
realized the significance of the Bernoulli’s work and used it as the basis, setting 
down the fundamental equations of motion and equation of countinuity for an 
ideal inviscid fluid in 1755, An alternate form of these equations was given by 
Lagrange in .1781 and 1789. Equations including viscosity were derived by 
Navier in 1822 and stokes in 1845. 

In 1858 Helmholtz published a paper on Vortex motion and in 
1 868 another on free stream line potential flow. 

Near the end of the nineteenth century. Fluid flow itself began to 
be extensively observed and investigated. The modern trend in fluid dynamics 
is to investigate the flow of electrically conducting fluid at very high 
temperature. 

At a microscopic level, the three states of a given substance are 
different because of the difference in the intermolecular distance. In many 
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There are two distinct methods of specifying the flow field 


cases at normal conditions the molecular distance of a fluid are less than the 
minutest of any physical dimensions of practical interest. As a result, we 
interested in the statistical average properties and the behavior of large 
numbers of the molecules, and not in that of individual molecules. That is, 
macroscopic and not microscopic, properties are of interest. In fluid dynamics 
as individual molecules are not being considered, the fluid can be regarded as 
a continuous medium and so the physical quantities such as mass, momentum 
etc. of the fluid contained in a very small volume are regarded as being spread 
uniformly throughout that volume. 


In dealing with the gases at very low pressure, as in the upper 
atmosphere, or at very high temperatures such as in a plasma. The continuum 
concept of fluid dynamics must be vilotated and the study that has to be based 
on the behavior of individual molecules (i.e. on the macroscopic approach) 


In our studies we frequently refer to a 'small element of fluid’ 
which is always supposed so large that it still contains a very large number of 
molecules, as fluid in a continuous medium. So when we take of an 
infinitesimal element of volume we always mean that which are physically 
infinite small. Such an element is called a fluid particle. 


in continuum dynamics we assume that the macroscopic fluid 
properties, for example mean density, mean pressure, very continuously with 
(a) the size of element of fluid considered (b) the position in the fluid, and (c) 
the time in (a), the variation becomes negligible as the element is physically 
very small. Thus, fluid properties density, pressure and velocity are expressed 
as continuous functions of position and time only. On this basis, it is possible to 
establish equation governing the motion of a fluid, which are independent of the 
nature of the particle structure. So gases and liquids may be treated together. 
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(i) LAGRAl^GIAM METHOD . 

in this method, the flow variables (velocity, pressure and density) of a 
selected fluid ^..lomient or particle are described. If To is the position of the center 
of mass of the fluid element at time to then the basic flow quantity in the 
Lagrangian description is the velocity V(ro, t). This method is also referred as 
‘individual time rate of change. 

(ii) EULERIAN METHOD 

In this method, the flow quantities are described at all points of space 
occupied by the fluid at all times i.e. flow quantities are defined as the function 
of position in space (r) and time (t). The basic flow quantity is the vector 
velocity V (r, t). This method corresponds to ‘local time rate of change’. 

Consider any scalar function f (x,y,z, t) associated with some property of 
the fluid (it could be density, velocity etc.) Suppose fluid particle has the 
position p (x,y,z) or p (r) at lime t. Keeping this point fixed the change in f is 
during the interval of time St is 

f(x, y, z, t + 5t) - f (x, y, z, t) 

Hence the 'local time rate of change' is given by 

(hf Lim f(x, y, z, t + St) - f (X, y, z, t) 

— = 5t->0 — - 

5t 5t 

As the point P is fixed the local time differential operator 5f / 6t is not 
carried along by the moving fluid. 

Now let at time t + 5t the fluid particle, which was at the position (x, y, z) 
originally Js in the position ( x + u6t), y + vSt, z + w5t) where u, v and w be the 
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velocity components at the position at time t. The corresponding change of f is 
given by. 

f(x + u5t , y -t- vat, z -I- wat, t +at) - f (X, y, z,t) 

or f(r+ar, t + at) - f(7, t) 
and rate of change is, 

df Urn f (r +5r, t +5t) - f (r,t) 

5t->0 

dt 5t 

This given tlie individual lime rale of change as the point p is moving t gives the 
rate of change which is carried along by the moving fluid. 

Since, 

f = f (x, y, z, t) 




df 

df 

df 

df 

df 

= 

— dx + 

— dy + 

dz + 

— -dt 
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9y 

dz 

at 
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af 

dx 
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dy 


af 

dz 


= 
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— . 


dt 


at 


ax 

dt 

ay 

dt 


az 

dt 


If V = [u,v,w] be the velocity of the fluid particle at p and 


dx dy dz 

... = U, -- = V, -.- = W 

dt dt dt 


Then 


df 

df 

df 

df 

df 

— = 

+ 

— u + 

V 

+ — w 

dt 

a 

dx 

dy 

dz 

df 

df 

df 

df 


+ 

( i -- + 

j- + 

k--) 

(ui + vj + 

at 

dx 

dy 

dz 


df 

df 




— = 

-- (V.V)f 




dt 

dt 





Where V = ui + vj + wk 

The first term on the right hand side represents the local rate of change of f and 
the second term the convective rate of change. 

If V is introduced for f in the above equation then the total derivative of 
velocity with respect to time is, 

dv dM 

a = -- = + V)V 

dt . 

In cartesion rectangular coordinates 

d d d 

(V. V) = u — + V — + w 

(3x dy Dz 

It is this term in the Euler relation for acceleration which is 
responsible for the non linearities in the equations of motion of fluid dynamics. 
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EQUATION OF STATE 


The electro dynamical state relation is a simple relation between 
the current density J, fields and fluid motion, or the ohm’s law is, 

J = a { E + pV + H } + q V 

Where a is electrical conductivity and q V is the current depending on the 
motion of the net charge q. 

The state of a compressible fluid is defined by the pressure P. the 
entsopy S, the internal energy E, the absolute temperature T. the mass density 
p and the specific volume V. 

An equation of state is formed by expressing only one of the 
thermodynamic variables in terms of the other two quantities mentioned above. 

for example, 

P = P ( P . T) 

or E = E(P, p) 

for a perfect fluid (no intermolecular interactions), having negligible 
viscosity the equation is 

P = RpT 

Where, R is the universal gas constant. Assuming the process to be adiabatic 
and isentropic we have, 

; ■ £ = CvT: : 
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= P / R p Cv as P = i> R T 


Also 

Cp — Cv “ R 
or 

Cp R 

Cv Cv 

or 

R 

Y-1=- 

Cv 

R 

Cv = — 
y-1 

P 

e = 

p(y-1) 

This is the caloric equation of state of the medium, is the internal energy per 
unit mass. 

When a fluid element changes its state isentropically, the state 

equation is 

d p 

---( --) = o 

dt 
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EQUATION OF CONTINUITY 


This ec|Liation simply expresses the law of conservation of mass. 

The quantity of fluid entering a certain volume in space must be 
balanced by that quantity leaving i.e. matter is neither created nor destroyed. 

We now formulate this principle mathematically. 

Let V be any arbitrary volume fixed in space, bounded by a 
surface S, and containing a fluid of density p. The volume element 6v is small 
so that p can be regarded as constant through it 



1 he mass of the liuid with in tiie volume V is Iv pdv the rate of generation of the 
fluid with in the volume is 
d 

— -Jv pdv 


dp 

J y — dv 

at 


for 5 5v 

p (dv) = pd(— -) = po 0 = 0 as volume is constant \A/ith respect 
Dt (It 

to time. 

if the volume V occupied by a moving fluid, the fluid enters V 
through parts of its boundary surface S and leaves through another part. 

A 

Let n be a unit outward normal vector drawn on the surface 
element ds. The normal velocity in n.v. The total outward flux is 

r , A - A - 

Js P ( p. V) ds = Js p. (() v) ds 

using Gauss's theorem 

= Jv V. (p v) dv 


The sum of the net outward convection of mass plus the rate of generation of 
the fluid with in llie volume must be zero. 

Op 

Jv V. (pV) dv + Jv — dv = 0 

at 

Op 

Jv[ — - + V. (pv)]dv = 0 

a 


or 


Since this is true for arbitrary elementary volumes 


dp 


a 


+ V. (p v) = 0 
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or 


f)C) d 

— _ + (p Vi) = 0 i = 1,2,3 

(3t OXj 


The Maxwell’s electromagnetic equations for a conducting medium are, 

DB 

Curl E = where B = |.i H 

a 

PD 

Curl H = 4 7t J + 

PI 

div B = 0 


and div D = 4 Tc q 


Differentiating the fourth Maxwell’s equation 

PD Pq 

div ( - — -) = 4 u - — 

a a 

using the second Maxwell’s equation it becomes, 

■ ' Pq 

div (Curl H -4 7t j )=4 TC — 

a 


Third Maxwell's equation is 


div H = 0 


we have 5q 

- 4 7T div J = 4Tt — 

at 


or 


5q __ 

- — -I- div J = 0 

at 

If the charged particle moves with velocity v 
J = qv then, 

— + div(qV) = 0 

at 

This is the equation of continuity for the electric charge moving under the effect 
of magnetic field which is similar to the equation. 

SHOCK WAVE, THEIR EXISTENCE 

If a small disturbance (i.e. with an infinitesimal amplitude and 
small velocity), is created within a non viscous isentropic compressible fluid it 
will propagate throughout the fluid as a wave motion and with tho velocity of 
sound relative to the fluid without suffering any distortion. 

In this case we get the wave equations, which shows that both 
the density and velocity variations follow the same wave patterns, derived from 
tlio lineal i/.iiiy ol the equations ol Iluid motion. I he solution (> (x,t) and u (x,t) 


both are the single valued function. The equations of fluid motion reduced to 
linear form as velocity is small, so that the term (v . grad) v in Euler’s equation 
of motion may be neglected. 

If the assumptions of infinitesimal amplitudes and gradients are 
removed, then the wave velocity will not be a constant and a simple wave 
velocity will distort as it propagates. 

In this case as velocity is not small, we have to solve the 
complete non linear equations of fluid motion. The method of characteristics 
can be applied to solve such hyperbolic type equations. 

A simple analytic solution of these equations of fluid motion may 
be find such that the density is a function of velocity only. After some 

manipulation we get equations in u and p. General solutions of these 
equations are, 

u = fi [x-(u±a)t] 

P = f2 [ X - ( u i a) t ] 

where fi and f2 are arbitrary functions, 
we consider the solution, 
u = fi [x-(u + a)t] 

This equation shows that the disturbance is propagated at an instantaneous 
velocity u + a, instantaneous because this velocity is a function of time. If the 
increment in velocity u is much smaller then speed of sound, then the solution 
IS the case of sound wave i.e. the linear case and the curve u = n (x - at) does 
not change its shape as the disturbance propagates, but if u is not small the 
shape of the wave is distorted as it propagates. 


The resulting distortions in velocity distribution are sho\A/n in 

figure. 



At a time ti, the points of high velocity move faster to the right than those of 
lower velocity. Thus the creast has moved faster to the right than the through 
and so the profile is disturbed. The compressive part of the wave, where the 
propagation velocity is a decreasing function of X, distorts to give a triple 
valued solution for u (x, t). which is physically impossible for longitudinal wave. 

Actually, the difficulty is over come by the formation of a shock 
wave. In such cases, a large variation in pressure and density occurs in a very 
narrow region in which the flow variables change rapidly and the fluid no longer 
undoiyoos isonliopic changes. Tho thickness, of this shock wave region is very 
small and so we may considered the shock wave as a surface of discontinuity 

for many piaclicnl problems of inviscid (liikis. 

The propagation of shock is faster than sound when observed 
from one side of the discontinuity and less than that of sound when observed 
from the other side. Hence the velocity of the shock will be supersonic viewed 
from ahead and it will be subsonic viewed from behind. 

The formation of these shock waves causes a great noise equal 
to that of an explosion, usually termed as supersonic bang so the effect of 
shock must be taken into account in the design of aeroplanes, pipe flow, 
supersonic flight of projectiles and so on. 


For the flow of an inviscicl and non conducting gas the laws of 
conservation of mass, momentum and energy are originally formulated in the 
differential equation form as it is assumed that the flow variables defining the 
flow are continuous functions. Flows are also possible, however, for which 
discontinuities in the distribution of these flow variables occur. So the 
conservation laws can also be applied to such discontinuous flow and hence 
across a shock. 

The condition for the existence of shock waves which may be 
called the jump conditions, relate the velocity pressure, density and 
temperature in front of the shock to those behind of them. The jump conditions 
are the simple consequence of the laws of conservation of mass, momentum 
and energy across the surface of discontinuity and the equation of state of the 
medium through which the shock is moving. 

SPHERICAL AND CYLINDRICAL SHOCK WAVE 

Consider the propagation of a shock wave, through a perfect gas, 
of great intensity (i.e. very strong) resulting from a strong explosion i.e. from the 
instantaneous release of a large quantity of energy. 

When the energy is suddenly released, in an infinitely 
concentrated form and distribution of density, pressure etc. depends only on 
the distance from some point then this is the case of spherical shock wave. 

When the energy is suddenly released along a line and 
distribution of all quantities in homogenous in same direction and has complete 
axial symmetry about that direction, then this is the case of cylindrical shock. 

Since we consider the symmetrical flow (centrally or axially) so 
there are only two independent variables, namely r and t. So the one - 
dimensional fundamental equations governing the adiabatic flow are 
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Where j = 0,1 or 2 stands for plane, cylindrical or spherical case, respectively, 
and r is the radial distance from the centre in spherical case, the radial distance 
from the line of explosion, in the cylindrical case and x coordinate in plane 
case; u is radial velocity: p the density. 

Similarity principle may be used to reduce these equations to ordinary 
differential equations. 

EQUATION OF MOTION AND JUMP CONDITION IN MAGNETO 
HYDRODYNAMICS 

In an electrically conducting fluids in the presence of magnetic 
field discontinuity i.e. shock wave in flow variable can exists. 
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The study of magneto hydrodynamic shock waves was begun in 
1 950 with the paper of F. de Hoffmann and Taylor. 

When electric current induced in the fluid, then their flow in the 
magnetic field produces mechanical forces which modify the motion, magneto 
gas dynamic owes its peculiar interest and difficulty to this interaction between 
the field and the motion. Thus the equations of magneto gas dynamic are the 
ordinary and electromagnetic equations. We only take into account the 
interaction between the motion and the magnetic field, we have ignored the 
Maxwell’s displacement currents. Through out the thesis we take p, = 1 
because magnetic permeability (p) differs only slightly from unity which is 
unimportant. We also assume that the dissipative mechanism such as 
viscosity, thermal conductivity and electrical resistant are absent. 

Since llio prohlcms (ionll in Ihis Ihosis rolalo to inngnolo gas 
dynamic or magneto - radiative shocks. We refer in this article the relevant flow 
and fluid equation. 



P = p RT ; 


Dp 

3U 

pu 



O 

II 

j 

1 

1 

I 

1 

Dt 

5r 

r 


Du 

1 

DP 


h 

ah 

ph' 

: + 

— 

— 

4- 

— 


+ — = 0 

Dt 

P 

dr 


P 

dr 

dr 


Dh 

Ph 

(lU 

Phu 


I- u— - 

-I- h + 

. = 0 

Dt 

dr 

dr 

T ■ . 
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d cl 

— (P p ') + u — (P p = 0 
Pt Pr 


where 

DP P 

Dt (It dr 

and u, p, p, h, r and t are the velocity, pressure, density, magnetic field 
transverse to the flow, radial distance and time respectively, T is the ratio of 
specific heat at constant volume and constant pressure, J=0, 1 and 2 
corresponding to plane, cylindrical and spherical respectively and u = 0 for 
plane and u = 1 corresponds to cylindrical and spherical both. 

In the presence of magnetic field, the relation connecting the flow 
and field quantities on the two sides of the shock surface are as follows, where 
the velocity in front the shock wave is zero 

h 2 (U-U 2 ) = hi u ; 

P 2 (U-U 2 ) = pi u ; 

1 1 
P 2 + — - h 2 ^ + p2 (u - U 2 )^ = Pi + --- hi^ + Pi u^; 

2 : ■ : . ■ 2 

1 rp 2 h 2 ^ 1 rpi hi^ 

— (u - U 2 )^ + — — + — - = — U^ + — + 

2 (r-1)p2 P 2 2 (r-1}pi pi 
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where suffix 1 and 2 correspond to the volume of the quantities just ahead and 
just behind the shock surface and u is the shock velocity. 

RADIATION PHENOMENON 

At very high temperature, the gases become ionized and radiation 
can be considered as a continuous emmision of energy is the form of 
electromagnetic wave which propagated in the medium with the speed of light. 
This energy is called radiant of thermal energy radiation where as according to 
ciuanluni theory the irrdiant onorgy omillorl or alrsorbod is not continuous, 
permitting all possible values, as demanded by the wave theory, but in a 
discrete quantified form, as integral multiples of an elementary quantum of 
energy, photon or light quanta. The amount of energy in each quanta being 
given by the product hu. where h is planck’s constant and u is the frequency of 
the radiation. 


Thus the quantum theory proposes the particles characteristic of 
radiation, while classical theory the wave characteristics, both being required to 
understand the complex behavior of radiation. A complete study of such a high 
temperature flow of a gas should consist of the study of gas dynamic field the 
electromagnetic field and thermal radiation flux simultaneously. The analysis of 
simultaneously effects and magnetic field forms the subject matter of radiation 
magneto hydrodynamics. This theory of thermal radiation can be applied to 
understand the processes which take place in stellar media, to explain the 
observed luminosity of stars and nuclear explosions and also to high 
temperature flow. 

Radiative transfer and radiative heat exchange have an influence 
on both the star and the motion of the fluid, this influence is caused by the fact 
that fluid loses or gains energy by emitting or observing heat. On the flow field 
of the gas there are three radiation effects expressed in terms of radiation 
pressure, radiation energy and radiation flux 
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RADIANT PRESSURE 


By the theory of electro dynamic the pressure of a radiation field 
is equal to one third of radiant energy is 


1 1 

Pr= --Er=-aRT"; 

3 3 

where T is temperature and an is Stefan - Boltzmann constant. This is the only 
component of the radiation which differs from zero. 


RADIATION ENERGY 

The radiant energy density Er per unit mass of the fluid is given 
by 

aRT" 

Er= 


where p is density of fluid. 

RADIATION FLUX 

The net amount of radiant energy passing through the surface per 
unit area per unit time is called the radiant flux through the surface is given by 

Cr OrT" 


or Fr = Dr grad Er 
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wlicic Dk == i-- lilt) l^osseland diffusion coofficicnl of radiation and L|< = — 

3 Kp 

the Rosseiand mean free path of radiation. 


Magnetic pressure Ph and the magnetic energy Ei, are given by 
h^ 

Ph = -- : 

2p 

h^ 

Eh= : 

2p 


where h is the magnetic field. 

Tlio fundamontal oqiiation of radialion inarjiiolo gas dynamics for ono 
dimensional flow are 




HI I) I I << 

-f p - — ( — ) -I- - ( Fr* ) = 0 

Dt Dt p pH d: 

where 

D a a 

— = ( — + u — ) 

Dt at ar 

where 

P = Em + Er + Eh 


Em = Pm/r -1 


and F is the radiation heat flux, where suffixes m, R, h attached to the 
quantities to expressions for material radiation and magnetic terms 
respectively. 

SiWILARITY PRINCIPLE AND SELF SIMILAR GAS MOTION 

It is not always possible to solve non-linear differential equations 
describing the physics of a motion or a process by using mathematical 
techniques. An approximate solution of such a problem can be obtained by 
solving a similar problem which is easier to solve. 

We consider the problem of one dimensional adiabatic flows of a 
perfect gas with constant specific heat, with either planer, cylindrical or 
spherical symmetry. The system of equation for flow of this type is given. These 
gas dynamic equations contain five dimensional quantities p, p, u, r and t. The 
dimensions of three of which are independent, for example density, distance 
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and time. So the equations admit three independent similarity transformation 
groups of quantities. Dimensional analysis can be used to obtain these 
groupings. By the successive application of these three groups we can obtain 
solution for the different flow similar to each other with altered density, length 
and time scales. Analogous transformations am made at the same time in the 
initial and boundary conditions of the problem. 

The motion itself may be described by the most general functions 
of the two variables r and t, p(r,t), p(r,t) and u (r,t). These functions also contain 
the parameters entering the initial and boundary conditions of the problem. 
They do not depend upon the position r and time t independently but are 
functions only of the combination (r / 1). In other words, the distribution of all 
quantities with respect to r change with time without changing their form they 
remain similar to themselves. This type of motion in which the distribution of 
the flow variables remain similar to themselves, (i.e, similarity in the motion 
itself) with time and very only as a result of changes in scale is called self 
similar. Consider the distribution of pressure. The function p (r.t.) can be written 
in the form p(r,t) = re (t) p (r/R), where n (t) is the scale of the pressure and R(t) 
is the length scale both depend on time in same manner, and the 
dimensionless ratio pin = p(r/R) is a function of new dimensionless coordinate r] 
= r/R. Multiplying the variables p and ^ the scalar functions 7i:(t) and R{t), we 
can obtain from the function p(ri), independent of time the true pressure 
distribution. The other flow variables, density and velocity, are expressed 
similarly. 

It is a natural question that what requirements must be satisfied 
by the conditions of a problem in order that the motion be self simitar. 

Dimensional analysis is used to answer this question. Since the 
dimensions of pressure and density contain the unit of mass, at least one of 
the parameters in the problem must also contain a unit of mass. In many cases 
this is the constant initial density of the gas po, which has the dimension ML'^. 
Let the parameter containing the unit of mass is a. It can be assumed that its 
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dimensions are [a] = mU'T®. The dimensions of the function p, p and u are [p] 

= ML'^ T [p] = M L'^, and [u] = LT we can without any loss of generality 

represent them in the form suggested by Sedov. 

a r 

— G, u = — V 

r t ^ t 

where p, G and V are dimensionless functions that depend on dimensionless 
groups containing r, t and the parameters of the problem. 

For self similar motion it is possible to reduce a system of partial 
differential equations to a system of ordinary differential equations for new 
reduced functions of the similarity variable q = r/R, R=R(t). The boundary and 
initial conditions of the problem are made dimensionless and in term 
transformed into conditions on the new unknown functions of -q. This simplifies 
the problem greatly from the mathematical stand point and in a number of 
cases makes it possible to find exact analytic solutions. 

The problem of a strong explosion represents a typical example 
of a self similar motion. This problem was formulated and solved by Sedov and 
succeeded in finding an exact analytic solution, to the equations of self similar 
motion. The same problem was also considered by Stanyukovich in his 
dissertation and by Taylor, both of whom formulated the equations for the 
problem and obtained numerical and not analytic solutions. 

The parameters in the problem of a strong explosion are the initial 
density of the gas po ~ ML'^ and the energy of explosion E ~ ML^T ■^. The 
energy E is always equal to the total energy of the moving gas, and as a result 
an energy integral appears in the problem. These two parameters can not be 
combined to yield scalar with dimensions at either length or time. Hence the 
motion will be self similar that is, will be function of a particular combination of 
the coordinate r(distance from the center of the explosion) and the time t. The 

23 ')' : 


-p. p = 


k 1 I 1 r, I ? 


m\(\a\ pressure and speed of sound p and c in llio problem of a strong explosion 
are assumed to be equal to zero, and hence these quantities are not 

[).ii iUiu ;1< M :i r>l tliis inobldin J'.o lh<! (|uaiitily i/l (.an not a(!iV(i aa tii(3 (iimilaiily 
variable. In tliis case Iho only dimonsional combination \A/hich contains only 
length and time is the ratio of E to po, with the dimension [E/po] = T Hence 
the dimensionless quantity 

Po 

n - r ( 


Can s(iiv(! as llu; s.imilaiity vaiialdo. Iho di;;liibution at pressure, density, and 
gas velocity can be expressed as functions of one dimensionless variable q. 

CONCEPT OF SELF GRAVITATION 


A fluid can be referred as self gravitating, when the fluid mass 
included is large and isolated and the gravitational attraction of other part of the 
fluid provides and volume force on any particular fluid element, for example, a 
gaseous star. 

In the case of spherical symmetry the effect of all masses on a 
particle (point mass) at a distance r from the center of symmetry is equal to the 
force of attraction by the point mass placed at the center of symmetry and 
having a mass m (r,t). 

The behavior of gravitationally interacting gaseous masses 
forming a star is given by the appropriate equations of motion. 

The continuity equation for radial gas motion with sphericaf 

symmetry is 
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for the inviscid gas, momentum equation taking gravitational forces into 
account (according to Newton’s Law of gravitation gravity g at r is g=Gm/r^) can 
be written as 


— + 


where G is the gravitational constant and m (r,t) is the mass of gas with in the 
sphere of radius r. 

To determine m we use the equation, 


4 7x p r^ 


which expresses the fact that the mass of a spherical shell of radius r and 
thickness 9r is 4 k r^ p dr and the energy equation for adiabatic gas motion. 

d p d p 

--( - ) + u --( - ) = 0 

5t p'*' 5r p^ 


These are the four equations with for unknowns p, p, V and m, 


A very large no. of papers have been published by now in which 
analogous self similar solution were obtained to explain the adiabatic unsteady 
flow in self gravitating gas and analysed for systems of partial differential 
equations encountered in various problems of astrophysics such as internal 
motion in stars, motion of nebulae etc. 
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CHAPTER - 2 


PROPAGATION OF PLANE RADIATIVE SHOCK WAVE IN 
UNSTABLE HOMOGENEOUS MEDIUM 

INTRODUCTION 


Scalar [2] and Pai [1] was the first person to investigate the 
problem of dynamical repture of the equilibrium of air masses in a non - 
radioactive state of a gas assuming that the velocity of an element of gas at a 
distance r at time t to be given by U=Gr/t where 0 is constant. Deb Ray [3] has 
studied the Same problem. Carrus Etal [5] discussed an analytical solution of 
same assuming that medium at rest. 


In this chapter we consider the disturbances headed by an 
radioactive shock surface of a rotating thunder in the earth medium and further 
we solve the problem of Deb-Ray [4], Further we assume that through out the 
flow the temperature gradient is zero and we take specific heat index y does not 
occur except that of term of internal heat energy. Which is a part of internal 
energy inside a thunder shock wave. In whole problem we take 0 = 2/3. As Deb 
Ray [4] in our problem also directly prove the result that mass enclosed within 
the surface of rotating shock wave at any time equals that under the same 
radius of a undisturbed state of the gaseous mass parcel of thunder system 
without imposing any special restriction on any of the parameters associated 
with the solution. The gaseous medium is taken to be in air. 

SELF SIMILAR FORMULATION 


The equation of motion governing the flow behind a shock wave 


are 



5p 

dp 

du 

_ + 

u— ~ + p ■ 


di 

'• dr 'v,':' 

/)r 


(2.1) 


Where u, p and p are the velocity density and pressure at the gas at a radial 
distance r from the origin at time t. m is the mass within a radius r and G the 
gravitational constant. T denotes the temperature. From equation (2.3) we get 

p/p = Constant (2.4) 

equation of continuity expressed as 


5m 5m 

— - + u =0 (2.5) 

5t 5r 


The disturbance is headed by an isothermal shock and the condition are 


P 1 (v-Ui) 

= po V 

(2.6) 

Pl-Po 

= poV Ui 

(2.7) 

Ti 

= To 

(2.8) 

mi 

= mo 

(2.9) 


Where po, po, mo and To are undisturbed pressure, density, mass and 
temperature just ahead of the shock and pi, pi and Ui are the quantities just 
behind the shock from equation 2.6 and 2.8 we get 



1 ^ 

and 

B =- (2.13) 

where = po/po is the square of the sound velocity ahead of the shock. 

Equation (2.6) shows that the mass particle behind the shock 
front moves with a velocity less than that of the shock front. Hence in the 
undisturbed state that mass confined with in the shock front at any time t must 
be equal to that contained within the surface of radius R. 

Where R is the shock radius at time t. 

i.e. R = R(t) 


'hus we have 


R R 

I 4 71 p r^dr =J 4 Tc poi^dr 
0 0 


(2.14) 


Equation (2.14) can be transformed to the nondimensional form 

1 p 

on, 

Where we have taken £= r / R 


and 


vr 


-n 


(2.15) 


V being a constant 
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Since the stellar model of infinite radius is initially in hydrostatic equilibrium due 
to its seif gravitation, 
we have 

Po = 2kGV^2-2^W(|i-1)(3-^) 

for the positively of the equilibrium pressure Po, we should have 

1<}x<3 (.216) 

let us seek a solution of equation (2.1) and (2.5) in the form 
r 

u=0 — (2.17) 

t 

SIlVlILARiTY SOLUTION 

In order to reduce the equation of motion to ordinary equation we 
now introduce the following similarity transformation, let similarity variables ^be 
in the form 


h 

= 

|.a 


h 


(n) 

(2.10) 

P 

= 

rk-.2pp(^^) 

(2.19) 

m 


r'^^^t^z(ri) 

(2.20) 


The constant K, X, a, and b are to be determined form the condition of the 
problem. 

The total energy Q inside a shock wave of radius R is given by 



R 1 

p Gmp 


Q -hr I (--pU^ 

1 )r^dr 

(2.21) 

0 2 

(y - 1) Y 


Where as the total energy with in the same volume in the equilibrium state 

given by 

R po 

GrnoP,, 


11== 4k 1 ( 

)r’clr 

(2.22) 

0 A 

Y-1 

Y 


Where 

R 

4 7t V r 


mo = 4Tcl por^dr = 


(2.23) 

0 

3 - jj. 



Since the motion is due to a straight excitation of an unstable condition of air 
mass we have 

Q = H (2.24) 

We consider the shock surface to the given ri by iio = constant 
Where Tjo is the shock front Value of ri 

We have V = - bR/at (2.25) 

Let us consider K=0, X,= -2 substituting equation (2.17)-(2.19) in the equation 
(2.10) - (2.12) and solving for Q and p (qo) and {r\o), we get 

b 

0 = - (2.26) 

a(1-p) 
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p[2 7xGv2/(|.i-1)(3-^)]R'^^^ 


V 

a(iio) = 

p R 

let p at the shock defined by 
Co^ ^ 

p = = R-^‘ r‘ 2 7x G V / (1.1 - 1) (3 - 

the right hand side of equation (2.27) — (2.29) are function of qo 

l-i a 

2 b 

thus we have without any loss of generality 

k=0, X- = -2, a = - a, b = 2 

hence we obtain from (2.26) 

2 

0 = 

h(1-P) 

which express the value of 0 in terms of p & p also from (2.30) 
we have 

P = Tc G V po / 2 (p - 1 ) (3 - p) 
which is a constant 

The following four ordinary differential equation are derived 
(2.1 -2.5) 


(2.28) 

(2.29) 

if 

(2.30) 

(2.31) 

(2.32) 

(2.33) 

from equation 



Inloini:'. ()fiil)vti-,in(|(?.17 ? ?1) mikI (P/U) 

- 0 + O^- -I- 1/L> (ti) [2 p (i]) - a n p (n)] + G Z{n) = o (2.34) 

p' (p)/pOl) = [4 + 21.1-5 |.L0/p (2 -p0)]1/ri (2.35) 

QSi/a(ii) = [(2-30)/(2-ia0)] 1/n (2.36) 

Z1 (p) / Z (p) = (2 - 30 ) / 2 - p 0 )] 1/p (2.37) 


After integrating (2.36 - 2.37) with respect to p we obtain the following result 


P = P,(,-/R)2-m(^-30/2-„0) 

(2.38) 

p = Pi (r/R) + <-' 

(2.39) 

m = mi(r/R) 

(2.40) 


where denotes the mass within the shock surface equation (39) we get 


J -- (2-p())/2(.1-p ) 

0 

I'l 

where as by (32) we get 
P / 3 - p = ( 2 - p 0 ) / 2 (3 - p ) 


(2.41) 


(2.42) 


So we find from equation (2.41) and (2.42) that equation (2.13) is identical 
satisfied for all constant values of 0, p and p subject to the shock condition 
(2.32) from (2.36) we also have. 


P(h) = P (ho) (p / ho) K^ + 2n-5a())/{2-^0)n] 

Q(p) = Q(po)(h/ho)'''''°^/'^-'‘°^’ 
Z(p)=Z(po)(p/po)'''-'">'''-^‘"” 


Using these values in (33) and (34) we get 
- 02 + o2 + 4 p (3 0 - 2) / p (2 - pO ) (p / po) '"“ + 
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(2.43) 

['.(jualion (2.^3) is to ho .also idonlically salisfy lor all v.aluo ol i) honco Iho Icmims 
involving i] either constant or zero. 

One set of solution is given by 

0 = 2/3 (2.44) 

from (33) by using the value of above value of 0 we get 
p = (3-^i)/3 (2.45) 

hence from (2.43) we have 

p - 48 + 36 = 0 ■ (2.46) 

We have on simplification 

p = 2.60 orp = 1.04 (2.47) 

both Ihoso v.aluos lies wilhin Iho tanrjf} given by ofiualion (2.1G) using abovo 
values of p the equation (2.45). 

We get 

[1 = 0.13 or |l = U.Gb (2.48) 

using the values for r, p, p and rn also for po Po and mo we get from (2.25) 

0" P / 4P(H -1) (3-n) [(2 - n9) / (10-2 n 0 -2n)l + (2 - |,i0) / (2p (y -1) (n -1) (3 - (i)") 

-2(2 - J.0) / (i(3 - M) (10 - 4n + pO) = [(2y - 1 - 2n (y-l)] / (y -1) (n - 1) (3 - n ) 
(5 - 2n) 

we have p, 5/2 and (j. 3, p. 1 (2.49) 
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CONCLUSION 


We find that the two possible sets of solutions given by 


0 = 2/3, p = 2.60, p = 0.13 and 
0 = 2/3, p= 1.04, p = 0.65 

where y has its independent value the equation satisfy by substitution in (2.49) 
the numerical values are different from the case of non radioactive besides 
these two sets there is no other possible solution of this two type in such 
system. Viz. from 0 = constant bearing any physical significance the value of 0 
i.e. 2/3 is the same for both the cases radioactive and non - radioactive. 

These two sets of solution are interesting one in which the pressure and 
density and uniform over that whole mass at any instant of time t (given). But in 
radioactive expansion cases the same fall in density as well as in pressure 
throughout the air mass of a sudden expansion in the atmosphere in the 
homogenous medium of gaseous is in an radioactive state as p / p = constant 
or p oc p ( p is directly proportional to pressure). 

These phenomenon are normally helpful in dealing with the 
astrophysical problems on account of clouds condition produce the plane 
shocks in the compressible gases of atmosphere. These are strong shocks 
which propagate with high mach number and dissipates in the lower strata of 
atmosphere on account of viscous layer prevailing. 
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CHAPTER -3 


ANALYTICAL SOLUTION OF SELF GRAVITATING GAS BEHIND 
A CYLINDRICAL SHOCK WAVE IN MAGNETO GAS DYNAMICS 

INTRODUCTION 

The unsteady motion of a large mass of gas followed by sudden 
release of energy result in novae and supernovae. The explanation and analysis 
of the internal motion in stars in one of the basic problems in astrophysics. 
Numerical solution for self similar adiabatic flow in self gravitating gas were 
obtained by Carrus etal at [1] and Sedov [3] starting with these results. Ryazanov 
[2] obtained a particular analytic solution. But this solution does not describe the 
flow behavior in general. Singh [4] has discussed the self similar adiabatic flow of 
self gravitating gas in ordinary gas dynamics and has obtained numerical 
solution. 

In this chapter we show that the azimuthal magnetic filed has a 
significant effect on physical parameters when the gas is self gravitating. The 
motion of the shock wave is assumed to satisfy the power low. 

po = AR-“ (3.1) 

where A and a are constants and R is the shock radius. The total energy of the 
flow increases with time because of the pressure exerted on the gas by an 
expanding surface. 

Therefore, 

E = B f (q > 0) (3.2) 

where B and q are constants and E is the total energy. The magnetic field 
disfribulion low is 


ho 


CR-l' 


(P>0) 


(3.3), 


wher© C and are constant and the value of q and p to be determined later. The 
magnetic field is directed tangential to the advancing shock front. 

The flow variable just ahead of the shock are 


4 Tt A 

Uo = 0, mo = R 


3- 


a 


and 


2 7tA2G 

Pq=: R2-2a ^ 

(a-1){3-a) 2p 


(3.4) 


where 


1 + (t = (X 


(3.5) 


and viscosity is neglected. 

We investigate three type of models, the first having total energy 
of the explosion to be constant (a = 2.5), the second having constant velocity of 
propagation of shock waves ( a - 2) and the third having neither constant total 
energy nor constant velocity propagation of shock wave (a = 1.5). An idealized 
magnetic field is considered for only a portion of sphere enclosing the origin i.e. 
the point of explosion. 

EQUATION OF MOTION AND BOUNDARY CONDITION 

The equation governing the cylindrically symmetrically flow of 
rotating gas under the influence of transverse magnetic field are written as. 
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wh6r6 r, t, m, u, p, p and y ara radial distanca from tha cantra, tima, mass 
contained in a sphere of radius y, velocity, density, pressure and ratio of two 
specific heats, respectively, as G represents the gravitational constant. 


The Rankine - Hungoniot boundary conditions at the shock (3.10) are 
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(3.14) 


111 


F 


mi 


me 


(3.15) 


where 


M' 


1 2(£-1) 1 y-1 

— + [ — + — ( e- 1) ] 

/i2 r^.., 4.'! \ _ r.\/ -"I ^ 1 K/i^ 4M^ 


yM 


[(y+i)-(y-i) ] M 


(3.16) 


dR 

V = — being the shock velocity and is give by the quadratic equation 

dt 


Ma^( 2- a) £^ + [(2y- 1)+ ]^-(Y + 1) = 0 


(3.17) 




The Mach number M and Alfven’s Mach number Ma are given by 






V^Po 

^2 _ gpj(j 

TPo 


2 _ 


V'po 

Ho" 


(3.18) 


let us seek the solution to the equation in the following form 
r 

u = - u (q) , p = r'^ t^ Q (q) 

t 


m 


r‘^'^^t^w(ri) 


(3.19) 


p = r t p (11), H = Ho = H. = r *'>«> ' " t ^ N (t,) 


(3.20) 


where iy = r^ t*^; while k, A,, a and b are constants and are to be determined from 
the condition of the problem. The total energy E inside the shock wave of radius 
R is given by 
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R 1 


E= 4n 1 ( -pu^ 


o 2 


p Grnp 

+ ... . cl r = B f ( q > 0 ) 

2 r 


y- '1 


(3.21) 


In terms of the variable q, we can express the total energy as 


471: Tio 1 ' p 

E = — I j-p (l<+5)/aH ^X.-2-(b/a)(k + 5) ^ y 2 Q .J. + ) 

a q 2 A-~1 




(3.22) 


Tlo and 11 being the values of q at the shock front ad expanding surface, 
respectively 

Equation (3.22) yields 


a -5 

b 4 + q 


(3.23) 


we choose qo to be constant at the shock surface. This choice fixes the velocity 
of the shock as. 


V = - 


b R 


a t 


(3.24) 


using equation (3.19) and (3.23) in equations (3.11) - (3.14) and in equation 
(3.18) we assume, without any loss of generality. 


10 


K = - a, A, = 0, b = 4 + q, a=-5anda 


4 + q 


(3.25) 
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Hence 


R= (326) 

where it is assumed that initially R = 0 from equation (3.26) it can be seen that 
the value q = 1 corresponds to uniform expansion at the sphere. 

ANALYTICAL SOLUTION OF THE PROBLEM 


Q’ 1 


5n 


( - 

- u) + (2 

- a)u 

- 5 q u’ 

= 0 


(3.27) 


LI 

a 








1 



P 


P' 

(/ - 3 

N' 

5u'ii ( 


■u) + 

u (u -1 ) 

+ — 

(1-(X)- 

- 5 q — 

( 

+ 5 q — ) + 


a 



Q 


O 

Q a 

N 


(2 

- ot) ( 

a-1) 1 





(3.28) 


n 


LW = 

0 







N' 1 

u - 1 


5q 

— ( 1 , 1 ) 4- LI - 5 q u’ + 

- (X u = 0 

(3.29) 


N (X 

1 



P’ 1 a 1 

Sq ( u) + u (1— (X ) — 1 - Y [ - (X u + 5 1] ' '-* ) 1 ~ 0 

"p (X LI a 


(3.30) 

w’ Q 

5q + 4 7t — -(2 - a) = 0 (3.31) 

"w "w 
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where 


and prime denotes the differentiation with respect to ii. The jump conditions at 
the shock are given by 


RESULT AND DISCUSSION 


Similarity solution of the problem of propagation of a spherical magneto 
gas dynamic shock wave have been obtained. For numerical calculation the 
flow and field variables have been taken in the following non dimensional 
forms. 
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CHAPTER- 4 



ANALYTICAL SOLUTION OF SPHERICAL SHOCK WAVE IN A 
ROTATING GAS WITH AXIAL COMPONENT OF MAGNETIC 
FIELD 


INTRODUCTION 


The propagation of hydromagnetic shock through a rotating gas being 
close to the actual situation of thunder. Kumar and Prakash [1 ,2,3] recently 
using C.C.W. [4,5,6] method have investigated the propagation of weak and 
strong diverging cylindrical hydromagnetic gas taking density distribution 
variable but axial component of magnetic field of constant strength. 


In this chapter, using the C.C.W. method we have considered the 
problems of the propagation of a converging and diverging spherical shock 
wave through a rotating gas under the influence of a magnetic field of constant 
axial, simultaneously for both weak and strong waves we have assumed initial 
distribution is constant. 


The analytical expression for shock velocity and shock strength 
have been obtained for weak shock under the two conditions namely when the 
magnetic field is weak and when the magnetic field is strong. For strong shock 
also we have considered two cases, i.e. when the magnetic field is strong and 
when po/pi is approximately equal to y+1/y-1 which is purely a non magnetic 
case. 


BASIC EQUATIONS BOUNDARY CONDITION AND AN ANALYTIC 
EXPRESSION FOR SHOCK VELOCITY 


The equation governing the spherical symmetrical flow of rotating 
gas under the influence of axial component of magnetic fields are written as 



where u, p and p are the velocity, pressure and density behind the shock front 

It'Mpi'clivi'iv W'hil*' 1 1. i‘i nvliil (MiiivitniK'itlM kI lii'I'l I I 


ut a sinyle puiamelei 


where 0 and 1 denote respectively the states immediately ahead and behind the 
shock front; U the shock velocity, a the sound speed (y Po/po) and the Alfven 

speed is (p Ho^ / po) 
where 

Ho^ = Hzo^ 

In the equilibrium state we have 
1 dpo 

po dr I 

Integrating the proceeding equation we get. 

A A 

P = Kilog - and ao = K 2 [log — f ^ 

t ' I 

where Ki = p’V^ 

yKi 

K2= (- — and 

P 

p’ log A is constant of integration. 

WEAK SHOCKS 

For a very weak shock we take the parameter N as 


N = H-£, 


£ < < 1 




Now coiisiilur Uic Iwo cuit.u'j ol wouk uiiU tilioiitj iiuiyuollc liold 


When the magnetic field is weak i.e 


and therefore, the boundary conditions [4.2] reduce to 


Case 2 : 

2 

When magnetic field is strong bo » 
circumstances the boundary condition (4.2) reduce to 


under these 
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TROf^G SHOCK 


When Ihe magnelic field is strong bo^ Qo' 
under t!iis condition tho boundary condition (4.2) becomes 


where 



For strong shock when bo « Oo ( i.e. when N -> 
mall the field may be regarded as independent of magnetic field 


For converging shock the characteristics form of system of 


equation (4.2) is 


where as the characteristic form for deverging shock is 


where c'' = a^ + b 


WEAK SHOCK WITH WEAK MAGNETIC FIELD 


Now siibstiluiing 



dpo dao 

^ , and 

Po ao YPo 

in equation (4.11) we get 

d£ [ y - 4] 1 A dr 

= r + log — ] - — (4.12) 

S 9 2 

Y Pi r r 

pHzo' 

where Pi^ = 

yki 

p/'-l) 

Integrating of equation (4.12) yields ■ 

1 A 

g =K 3 r[T- 41 /yl^^ 2 exp{ (log -)'} (4.13a) 

2 (^ 1 ^ r 

where log K 3 is constant of integration in case of diverging weak shock with 

weak magnetic field, we have the relation 

A 3 A 

£ = K 3 r ( log — ) exp { - Pi^( log ) — } (4.13b) 

r 8 r 
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where log K 3 is constant of integration and 


3 

qi = ( — + — ) 

4 4 

Hence from equation (1 3a) and (1 3b) we have 

u Y ^ ^ 

... = 1 + K 3 r''-<'-'V>-‘“''Vexp{ (log -f (4.14) 

ao 4 2P,^ r 

y y +1 A 3 ^ 

= 1 + K 3 r ( log — - )''’i exp{ - -- (log — )’^ (4.15) 

ao 4 r 8 r 

for converging and diverging shock respectively. 

WEAK SHOCK WITH STRONG MAGNETIC FIELD 

Substituting the shock conditions (4.7) in (4.9) we get 

dpo pHo^ dr dbo 2 dpo 

d£ +{ — - (--- + — )“ } =0 

Po YPo r bo YPo 

Substituting the respective quantities in equation (4,16) we get 

d^ y(4+Pi^-2) a dr 

log — )■■' — (4-17) 

i 

y ' f 
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integrating equation (4.17) we get 


where log K 4 is constant of integration and 


when we consider diverging shock with strong magnetic field we get 


where log K^t-. is constant of integration equation (4.18a) and (4.18b) disclose 


for converting and driverging shocks respectively 


STRONG SHOCK 


Substituting shock condition (4.8) into (4.9) we get 


4 “ 


( -r + U^[- 

y 2 N Y po y 


dpo yki f3i^ (N-1) + VNx dr N VNx 

Po Po (N-1)-VNx r [(N-1)- 


Substituting the respective values we get. 


du^ B 


+ .„ u^-cr'^ 


dr r 


where B = BVM 


C = CVM 


X(N-1) 


B’ = 


(N-1WXN 


X (N-1) X 
M= 

y 2 N 


and 


Lki ykiN^ (N-1)+WNX 

C’ = — — + { ——————— } 

p' p’ (N-1) WNX 


0 (N-1)WNx 


} = 0 (4.21) 

VNx] 


(4.22 


N^NX 


(N -1) - VNX 




respectively 


In the last, the expressions for the velocity, the density, and the pressure of the 
gas J ust behind the shock surface are 


Integration of (4.22) yields 

C 


= Ksr'* + 


(4.23) 


where K5 is constant of integration. In the case of diverging strong shock with 
strong niagnolic field wo have 


U2 = Ksr-i' + 


(4.24) 


Where K5 is constant of integration finally, for converging and diverging shock 
we have 


(K5r;®+ -- ) 
B 


(4.25) 


ao' K2' log -- 
r 


and 


(4.26) 


WEAK SHOCK WITH WEAK MAGNETIC FIELD 


A 


A 


U = K2K3 ( log r ^ I^g ^ )2 ^ (4.27a) 

r 2Pi2 r 


ASA 
U = K2K3 log exp { - — log ( — )'^ } 

r 8 r 


(4,27b) 


A 


p= p’ [ 1 + K3 exp {- -——(log —f (4.28a) 

2Pi2 r 


A 3 A 

2/ 1^^ N-1 


p = p’ [ 1 + K3 (log — )'Sexp { - — pi^( log — ^ (4.28b) 

r 8 r 


A 


A 


p = ki(log - >Yk3r^^-"^'^V(log -)exp{ - 


2pi' 


A 

(log -f) 


(4.29a) 


A A 3 A 

k-i ( log — + y ks (log — )^ '’1 exp {- --pv (log — ) } 
r r 8 r 


(4.29b) 
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WEAK SHOCK WAVE Wl lil S I KONG WlAGNEIIC ! lELD 


A 


U= piK2K4(log -r2 

r 


A A 

2 o 2 \ \“1 


U = P 1 K 2 K 4 ( log - exp {- 2 Pi 2 (1 - Pi") (log -r } 

r r 


A 

p = p’ [ "I K4(log -- ) ^2 ] 

r 


(^= p 


A A 

[ 1 + K 4 (log -- exp {- 2 Pi^ { 1 -lK^) log ( — } 


A A 

p •“ l\| I I' Ml ’ V 'lAi ( l<ni )' I 

r r 


A 


A 


A 


p = kt [ log — M- Y K 4 ( log — )''i^ ] exp { - 2 (Vi { 1 - (K ') log(' 


STRONG SHOCK WITH STRONG MAGNETIC FIELD 

(N-1) C 

U= --- [K5r'''+ -}^ 


N 


B 


(4.30a) 

(4.30b) 

(4.31a) 

(4.31b) 

M n?i») 

)'M] 

(4.32b) 

(4.33a) 
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N for converging diverging cases 


where section (a) and (b) denote the expression for converging and diverging 
shocks respectively. 
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CHAPTER - 5 


ANALYSIS OF MAGNETO HYDRO DYNAMICS SPHERICALLY 
SHOCK WAVE IN A MEDIUM WITH VARIABLE DENSITY 

INTRODUCTION 

Pai [1], Green Span [2] Chester and Helliwell [3] have studied a 
cylindrical shock wave of line explosion, produced on account of a sudden 
release of a finite amount of energy expanding outwards in conducting gas 
subject to magnetic field. Deb Ray [5] and Singh & Vishwakarma [6] have 
obtained a similarity solution for strong cylindrical shock wave in a conducting 
non uniform medium Singh & Srivastava [6], have obtained the similarity solution 
for a axially symmetric explosion model in magneto gas dynamics I & II without 
thermal radiation. Singh & Pandey [7] discussed the conduction for the total 
energy of the cylindrical shock wave to be assumed as constant. 

The present chapter deals with the problem in which the amount of 
energy is not finite in an infinite concentrated form is suddenly released with 
spherical symmetry in a conducting gas subjected to an azimuthal magnetic field 
which produce a spherical shock wave. This propagates in the medium whose 
density varies with time. The motion of the gas is assumed to be isothermal and 
total energy of the wave. 

EQUATION OF MOTION AND BOUNDARY CONDITION 

The equation of motion for one dimensional unsteady isothermal 
flow of a perfect gas with transverse magnetic field are 

dp dp du u 

+ u + p + 2p — = 0 (5.1) 

dt dr Dr r 


au 

du 1 

dp 

+ 

Li + - 

H 

at 

Pr p 

hr 

r)h 

ah 

Pu 

-f- 

II + 

h — - 

f3t 

Sr 

dr 

aT 



— = 0 



Pr 




.. = 0 
p Dr pr 


(5.2) 


u 

-I- h — = 0 
r 


(5.3) 


(5.4) 


The magnetic permeability is taken to be unity other symbol have their usual 
meanings and they are function of r and t only. Where V is ratio of specific heat. 


Following Whitham(1959), we find the shock boundary conditions 


to be 


1 

u,- (1 ) 

N 


pi ~ Npo 
Hi = Nho 

P “ Pl^PoPo 

+ (2d + 1) N-2 Ma^ = 0 
1 Ma 

D - - ( -- )^ 

Y M 


(5.5) 


(5.6) 

(5.7) 

(5.8) 
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mimbo!. 

SIMILARITY SOLUTION 


Let the similarity variable ii be in the form 
r 

n= R = R(t) 

R 


(5.9) 


We seek the solution of the partial differential equation (1 - u) in the form 


u = RV in) , 


p = po g (n) 


p = po f(ii) 


h = R H (ri) 


(5.10) 


where v, g, f and H are function of r| the scale R, po and R time dependent in 
some manner as yet unknown let us assumed the radius and density of the 
cylindrical shock front are 


R = A t" 


(5.11) 


po 




(5.12) 


we now substilutG the relation (5.10) in to continuous (5.1 - 5.4) taking account 
of the definition of similarity variable (5.9) and applying relations of the type. 


9p R 

= gpo' - pog’n — ■ 

a R 


(5.13) 
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Oh H 2 

....= (2 R po + R Po) - Po''^ R’ 

D[ 2po Fl'lH’ 

Op Po g’ 

_... = — ( 5 - 15 ) 

Or R 

5h R’ 

_ „ 1/2 (5.16) 

- po 

OH RH’ 

to transform the derivatives the dot represent the differentiation of scales with 
respect to time and prime denote differentiation of functions with respect to 
similarity variable therefore the non differential form of ordinary differential 

equation and are 

(3 v’+(v-ii)g’ 2v 

+ - = 0 ( 5 - 17 ) 

a g h 

2(a-1) + {3 (v-h) H’ 

...... + V = 0 (5-1 8) 

2a H 

1 1 H 

(X - 1 ‘ 

-).v + (v-i])V'+ — {f' + H(H‘ + — )} (8.19) 

cc g 
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f' g’ 

= 0 

f g 

consider now the total energy E of the flow as 

R 1 p h’" 

E = 2nl [ + — -+ --- ]r.dr 

° 2pv^ y-l 2 

Equation (5.10) reduce the above equation in the form. 

11 f H2 

E = 2 n p t \ ( + + — ) ■n 

i-lo 2gv^ y-^ 2 

The total energy of the wave is constant so we have 
p = 2--4(x 


( 5 . 20 ) 


(5.21) 


(5.22) 


(5.23) 


it is clear from (5.23) that uniform expansion of wave take place, when p = 2 the 

solution have their physical significance only 

when 

-2< = p< = 0 

after using similarity transformation (5.10) in the jump condition (5.5 - 5.8) at the 
shock from yield 

g(1) = f(1) = N 


H(1) 


N 


Ma 


(5.24) 


y (1)T= (1 -1/N) 
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UBiiiij !<.K.(3‘s iiialliud btaiiiiuj Iioiu slujck lionl wlUi ua 
condition and counting up to time explosion equation are i 
calculated our result in the form. 


ADIABATIC CASE 


as usual the magnetic permeability of the medium is taken to be unity following 
with an (5.1) the shock boundary condition are taken to be 


Ill v\i h 


hi = \i/ h 


1 2(vi<--y) 1 (Y-1)(r'1) 


p,, = t ..... + -"--}lPoF^ 

(Y+1)(y-1) 4 Ma^ 


i2 


(5.27) 


where v can be calculated from quardratic equation 


2-y 2 Y 

[ ] i|, + [ — + + (Y + 1) ] \|/ + ( Y + 1) = 0 

Ma' 


(5.28) 


2 Ma^ 


and 






Ma- - 


M' 


(5.29) 


ir^ 


YPo 


SIMILARITY SOLUTION 


It is convenient to assume shock radius and density varying as time 


R = C e"^‘ , p = D e 


nt 


(5.30) 


Where m, n are constants. 


Let the similarity variable ii be in the form 
r 

,, = R = R(l) 


R 
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We seek this solution of partial different equations (5.1 to 5.3) and (5.26) in the 
form given in (5.10) 


where v, g, f and H are function of ri 


We now substitute relation (5.10) in the equation (5.1 - 5.3) and (5.26) and 
taking up to account of similarity variable (5.31) and applying the relation 
(5.31 - 5.16) we find the dimensional form of ordinary differential equation which 
yield by considering the possibility is 


Constant 


Constant 


The time dependent equations (on account of parameteni) are 



where dash denotes differentiation with respect to “ii” considering the total 
energy E of the flow as 


Equation 5.30, 5.31 and 5.10 reduce the above equation in the form 


In order to have total energy of the wave is constant. We have 


It is observe that uniform expansion of wave take place 


When 


riio soluUon their pliysicul significance only when 


-4 < = n < = 0 (5.39) 

After using similarity transformation jump conditions at shock proof yield. 
g(1) = vy 

v(i) = (1 _1/v|/) (5.40) 

f(1) = kv|; 

H(1) = m/ZMa 

Using R.K.G, Method from shock front with transformed boundary conditions and 
containing up to line of explosion for a adiabatic case equation (5.32) — (5.35) 

are integrated numerically on R. K. G’s program in the 




RESULT AND DISCUSSION 


The Equation of motion from (5.1 — 5.4) under the given boundary 
condition (5.5 - 5.8) are integrated with the help of (5.13 - 5.16) which under the 
similarity condition given equation number (5.10) which yield to non dimensional 
ordinary differential equation (5.7 - 5.20) for the propagation of Magneto Hydro 
Dynamics spherical shock wave. In an unsteady isothermal flow behind the 

perfect gas. 


The numerical values are plotted against the graph by fig ($•') to fig 
(?'M) for the various values of N = 10 , 100 the graph clearly indicates that the 
pressure, velocity, density distribution varies exponentially as we move away 
from the shock surface where as magnetic distribution tends for N = 10 it limits 
indefinite range 0.1 to 0.02 again it exponentially for the large value of N = 100. 
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CHAPTER -6 


PROPAGATION OF PLANE SHOCK WAVE IN MAGNETO GAS 

DYNAMICS 

INTRODUCTION 

The phenomena associated with heat transfer in rotating fluids 
are extremely complex whether the fluid is at rest or in motion, study or unstudy 
Wang [9] considered the piston problem with thermal radiation for one 
dimensional unstudy shocks using the similar method of Sedov [6], Hellwel [2] 
look a more general case of piston problem with radiation heat transfer for 
general optically and transparent limit. Elliot [1] discussed the self similar 
solution for spherical blast waves in air using the Rossland diffusion 
approximation under the assumption that there is no effect of heat flux at the 
centre of symmetry laumnch and Probrteim [4] have studied the radiation 
effects on shock waves in an spherical medium several other author e.g. 
summer [8] magnetic field also be taken into account. The unsteady model of 
Roche - Cosits of a gas distributed with spherical symmetry around a nucleous 
of large mass. It is assumed that the gravitating effect of the gas neglected 
compared with the attraction of the heavy nucleus Rosenau [5]. 

In this chapter we consider propagation of plane shock wave 
under isothermal condition. In the isothermal condition temperature gradient 
becomes zero behind the shock and radiation effect are already implicity 
involved. An idealized magnetic field is considered for only a portion of sphere 
including the origin i.e. the point of explosion Summer [8] the field in directly 
tangential to the advancing shock front. The shock is assumed to advance into 
conducting gas of specially decreasing density and produced by a specially 
decreasing magnetic field the total energy of the above in non - constant. 


equation of motion for adiabatic flow 



a,. 

au 


+• 

U 

+ p — 

= 0 

fTt 

ar 



8u 

au 1 

ap h 

ah h^ 

--- ■'I* u 

— + - 

+ — 

+ — - 

5t 

ar p 

ar p 

ar ar 

9P 

0p 

ap 

af 

— _ + u 

— + Y p 

— (y-1) 

— = 0 

a 

ar 

ar 

ar 

ah 

ah 

au u 


+ LI 

— + h 

— + h — 

^ r: c 

at 

ar 

ar r 



+ = 0 


(6.1) 


(6.2) 


(6.3) 


(6.4) 


Where u, p, p, h and f are the velocity, pressure, dens.ty, 

.diative heat flux at a radia, distance r tron, the center of the core - 1 - 
denotes the constant mass of the core, G be the gravitational constant, 

magnetic permeability Of the medium has taken to emity 

The equation of state for ideal gas is given by 


P = rpT 


(6.5) 


Where F is the gas constant. 

Also assuming local thermodvnantios eguiiibrium and taKen Rossland's 

diffusion approximation, we have 
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m 


a 


- {aV) 


( 6 . 6 ) 


3 ar 


Where (1/4) ac is the Stefan - Boitz man constant, C is the velocity of light and 
the mean free path of radiation, is a function of density and absolute 


temperature T, following Wang K. C. (9) 


n = po p 


‘j\\ 


(6.7) 


where no, a, and p being constant, the total energy of the wave is dependent 
on time as 


E = Eo f 


( 0 < q < 4/3) 


( 6 . 8 ) 


Where Eo being the constant in the front of shock wave, we assume that 


Pi 


= AR 


( 0 < (0 < 2) 


(6.9) 


where A is constant and R denotes the radius of the shock surface ahead the 

shock, the magnetic field distribution is assumed to be as 


hi = CR 


-k 


(2K = 03 + 1) 


( 6 . 10 ) 


At the equilibrium state the other flow variable ahead the shock are 


Ur-0, pi = [ 


AmG 

^ R -(1+«) 4. — (6.11) 

1+0) Kfl-k). R'* 
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the flow variable just behind the shock denotes the mass 
K) shock and v bo Iho volocily of Ihe shock and is given 


Where suffix 2 denote 
lliix per unit, across II 


SOLUTION OF THE EQUATION OF MOTION 


where C. W and K are constant the disturbance is heated by an isothermal 
shock and Ranking Hugoniot jump condition are. 


(G,12) 




Where V is constant to be determined by the condition that assume the value 1 
the shock front hence we have 


By use of equation (6.18) in equation (6.6) with the add of equation (6.5) we 
obtain 


Equation (6.1 - 6.4) and equation 6.6 are then transformed with the relation 
( 6 . 18 ) and (6.20) to the form. 



where L = ( 


is a dimensionless radiation parameter 


EQUATION OF MOTION FOR ISOTHERMAL FLOW 


The equation of continuity, momentum and magnetic field are the same 
as in the adiabatic flow but the energy equation becomes 


T = Constant 


By use of equation[4-^we have 
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For the isothermal flow we also assume that the jump condition are heated by 

isolhormni flow ns given by oquntion 6.12 - 6.10. 


Ify niihsliliilioii of lolalion (0,18), (0.20) nnd (0.22) in ocitinlion (0.1 - 0.3) ns 
(0.28) llin followinf) not of oqnntion nro nlilninod 


(6.30) 


RESULT AND DISCUSSION 


The transformed jump condition at the shock for both adiabatic and isothermal 
flow are given by 

2 

V(l) - (-)ii (0.35) 

3 

V 

D (1) = (6.36) 

(1-Q) 

4v a + 1 Q (2 - Q)M^a 

P(i) = — [ + ] (6.37) 

9 yMV {2 ( 1 - 0)2} 

2 1 

N(1) = -Vv( --)Ma (6.38) 

3 1-0 

and 

8v (1-0 )yM 2-1 M2A(yO-2)0 (2-0)0 

Q(1) = — . yQ { } + { — — — + — — ] (6.39) 

37 (y-1)M2 2(1-0)(y-1) 2 

Where V is obtained from the relation 
9 m2 

V= ( ---) v M2 [2 ( 1 +0)-y(-'“-) -n)]'^ (6.40) 

2 Ma 
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For oxhibiliiK] Hio niimorical nnhitinn, il in rniwnninnt In wiito vnriablo in nnn 
cliineiinional fuiiii an. 


u v(ii) 

--- = F[ -] (6.41) 

U2 V(1) 

P 1 D(ti) 

- = i ] (6.42) 

pj V’*’ D(1) 

P I P(ll) 

-= -- [---] (6.43) 

P2 if P(1) 

h 1 N (ii) 

...= (6.44) 

h, 11(1""'' N (1) 

and 

F 1 Q(ii) 

[„-] (6.45) 

F2 rf Q(1) 

The numerical integration is carried out for boath adiabatic and isothermal flow 
separately on DEC - system 1090 computer installed at IIT Kanpur by the well 
known RKG system programmer for the there values W = 1.00, W = 1.50 and 
W = 2.00. The other constants are = 20, M^a = 10, 6 = 2/3, and L = 10. 


The nature of flow and field variable for both the adiabatic and 
isothermal cases are illustrated through figured -^Hit is clear from the figure 
that velocity, density, pressure and magnetic field are maximum at the shock 
front and decreases rapidly towards the centre of explosion, when the flow 
behind the shock wave is isothermal but these variables are minimum at the 
shock front and increases steadly so we move towards the centre of explosion 
in adiabatic case. The condition of radiation heat flux has been illustrated in 
figure 5 we therefore conclude that the flow and field parameters are being 
subsidied, when the assumption of adiabatic to be not valid and the 
temperature gradient of the flow becomes zero. 





Without magnetic field 
With magnetic field 


Velocity Distribution 
Fig (6.2) 
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CHAPTER -7 


RADIATIVE EFFECT OF POINT EXPLOSION IN STELLAR 
BODIES 

INTRODUCTION 

Cnirus olnl [1] nncl Sndov |2] worn first to discuse tho bodies of stellar 
oxplosion in whicli n star nonnidorod In bo a potfnci self gravifoting gaa. 

In this chapter we consider the some point explosion in stellar bodies 
since the temperature of material is very high the radiation effect can not be 
ignore. The aim In this cliapter is to discuse the equation for the propagation of 
a radiative blast wave produced by a sudden point explosion in self gravitating 
system such as starts. We have taken similarity solution and medium has been 
assumed to be perfectly conducting plasma with radiative parameters 
independent of magnetic field. The disturbance is bounded on the outside. In 
this case also the viscosity and heat conduction are neglected and it is 
assumed that the flow is isentropic along a stream line, Lee [4] discussed the 
model of blast wave model to account for interior energy. In this chapter we 
also consider effect of transvers magnetic field in stellar bodies. 

EQUATION OF MOTION 

The equation of motion, continuity, energy and the filed equation 
in the case of a radiative gas are 

9u 5u 1 8p Ho d Ho Gm 

— . + u ___ + — — + — + — — s 0 (7.1) 

M dr p dr p r r^ 


du 2u 

dl I 


dp 

dl 


dp 


I, 


{7.2) 



5E 

5E 

0 

0 1 0(fr2) 


+ 

u — + 

P(- 

•(1/p)+u — (1/p) + = 0 

(7.3) 

5t 

5r 

r7t 

0r pr^ 0r 


5Hi) 

5H() 


0u 2u 


■I- 

U i 

llo( 

- ' - ) - 0 

(7.4) 

5t 

5r 


0r r 


where 

Ho" 


Ho" 


E 

2p 

P = 

2 

(7.5) 


The quantities u. p and p are radial velocity, pressure and density at a distance 
r at any time from the point of explosion the magnetic field has component (0, 
0,Ho), f is the radiation flux. G is the gravitational constant and m is the mass 
within the shock front at any time t such that 


(7.6) 


(7.7) 


(7.8) 


where i is the cofficient of opacity, C the velocity of light, 
we assume as in [5] 

P =ZP, (0<Z< 1) 


5m 

5r 


47Tp r 


We have 


E = 


and 


p(r-i) 

p dr 
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So that 


P 

E = (7.9) 

p(k-i) 

where K is called Klimishin’s cofficient and is given by [5] 

4(y-1) + Z {4-3y) 

K= (7.10) 

3(y-1) + Z(G-3y) 

Y being the usual ratio of specific heats, with the help of (7.9) , (7.2) and (7.4) 
and (7.3) can be written as 

ap dp au 2u k-1 a(fr^) 

+ u— + Kp( — + — )+ =0 (7.11) 

at ar ar r ^ dr 

let the motion be assumed to be confined within the shock front at r = R(t) then, 
the velocity of shock moving outwards is given by 
dR 

V= (7.12) 

dt 

SELF SIMILAR VARIABLES 

The following similarity forms are used for flow and field variables. 

P= PoR'^f(h) (7.16) 

p= po M/ (n) (7-14) 

u- R" (p (n) (7-15) 

H= Ho R'^g (h) (7 16) 

F = fo R' n (n) (7 17) 

m = mo R m (n) ^ 
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^I^ere ti=i 7R is a nondimensional radial variable f, \\f, (p, g, m, n are function of n 
only and tlie suffix zero denotes quantities in the undisturbed state, putting the 

equation (7.1) - (7.4) and (7.6) into their similarity fornn, we get 

R / R'^ (ricp - ri(p’) + 9’ 9 = 1/po9 


p^R(m-2n)F’ + HoR' 2 ^“ 2 "^g’ + moR 


(b-2n-1) 


m ( Po\p) 






R 

(n(p’+ 29) = P ( n 9) 

K'^ 


R 

R" 


29 f (k-l)fo 

h 


2n 


(mf - iVf2) = K( 9’f + - ) + ---- R (n n’ + -- ) 


l-m -n 

“V “ 

Po h 


(7.19) 

(7.20) 


(7.21) 


R g’ g 

— (r,k - Tj^ — ) + 1191 + h9’ + 29 = 0 

R" g g 


m 


4tcPc 


mo 




(7.22) 


(7.23) 


in order that all the unknowns may be expressible as function of ri along the 
following relation must be fulfilled. 


m = 2n, k=n, L=3n, 

and 

R 

= C 

R" 


b=2n+1 


(7.24) 


(7.25) 


Integrating (7 25) we get 



1 being an arbatrary parameter and C an absolute constant, As t->0 and R->0 
ukI tio w« imuil havt! n- I aiul A-0, lt\<}ii tilKXik ladlun l< Ih (jlvtsii as 


with the lielp (7.24) (7.25) and (7.26) the equation (7.19) - (7.23) become 


2(pf (k-1)f( 
C (2ii f - nf ) = K (f(p’ + ) + 


P = Po (^l) 


H = Ho R"g (ii) 


INITIAL CONDITION AND SOLUTIONS 




Let the explosion take place at a point at t=0 then, 


u 

zz 

CO 

for 

r 0 

u 


0 

for 

r--A 0 


and also as obtained earlier R->0 for any admissible solution, p, p, H and u are 
all finite for t— > 0 everywhere in r < R by putting the value of n=-3/2 in (7.26) we 
get the shock radius R as 

R = p P (7.32) 


p is an absolute constant from (7.32) the shock velocity is then given by 


dR 

V- — 
dt 

2 

= --(R/t) 

5 

(7.33) 

Again by using (7.32) the similarity transformation can be written as 



u = t-''®f2(J;) 


P= t-®'®f3(^) 

H = t U © 


F= f®^®f5(^) 

m = t fe © 

(7.34) 

rt-2^^ 


(7.35) 

The energy equation can be written in the form 


dE 1 

m r^ 

d 

— (UIt + F) = 0 

dr 

(7.36) 

E = %pu^ + 

P Gmp 

(K-1) 2 r 

(7.37) 
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Kp 


Gmp 


lr = y 2 pLi^+ + Ho (7.38) 

(k-1) r 

puUiog (7.37) into Ihoir similarity form we get 
2 

E = y2f(^)f®^®f2(^) + (1/k-1)f®^® 

F3(^) + y2t-®'®f4(^) (7.39) 

or 

E = f (0 
where 

1 f6(^) f(0 

F(t;) = y2 fi f2(^) + f3 (ty) + ya f4(^i) (7.40) 

k-1 (^) 

from (7.39) we have 


aE 


■ =f(^)f®^^ 


(7.41) 

ar 

aE 

6 

2 



. = „ 

t-ii'5f (^)_t-6/5f (^)— rf^^® 

(7.42) 

at 

5 

5 


with the help of (7.39), (7.41) , (7.42) we get 


aE 

6 

2r aE 



- + —Et ' + 

= 0 

(7.43) 

at 

5t 

5t ar 


aE 

2 

ar^ 




— E =0 

(7.44) 

at 

5tr^ 

-'dr ' 
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from (7.44) we have 


d 2rE 


+ (UIt+ F) = 


(7,45) 


5r 5t 


(U1t+ F) 


where we have taken the constant of integration to be zero and the value of V 
is substituted from (7.33). The equation (7.45) can also be written as 


x(^^^ pu^ + 


Gmp 


2 Rx 


t1/2 pu^^ + 


Gmp F 


where U = vu’ , 
from (7.46) we have 


r= Rx 


(7.46) 


p k-1 Cam 

= [ciu^^ (x-u’) + Cap (x-2u’) + (u’-x) - C4 f] (7.47) 

p (ku'-x) X 

where ci, Ca, C3 and C4 are function of time and H = cp as a consquence of (7.2) 
and (7.3) from (7.2) and (7.1 1) we obtain 


1 ap 

p dr 


(k-1) Sp 1 ap (k-1) 3p k-1 5(fr^) 1 au 2 


p ar pu at au at r ar 


u ar r 


(7.48) 


by putting the value of ^ and v in (7.34) we get the relations 


ap rv 3] 


(7.49) 


at R ar 
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integrating this we get 


substituting (7.47) in (7.53) we have 


Eliminating P between (7.53) and (7.47) and dropping the primes 
We obtain the equation 


which determines p. For simplicity we can write total derivatives in place at 
partial ones. Di, D2, D3, D4 and D5 being constants depending upon time 


flp 3pv rv dp 

at R ar dr 

by using (7.49) and (7.51) we can write equation (7.48) 


(7.50) 


1 5u 1 

V ar R (k-1)u 

u r u r 

( ) v( ) r^ 

V R V R 


afr^ 

ar 

(7.51) 


Similarly eliminating p between (7.53) and (7.47) we can get 


(k-2)/ Ai (y-2u) 2/k-1 


Az (k-1)ue d(fy^) 

(ku-Y)-{ I dx} 

Y^(u-y) ""p (u-y)y^ dy 


(A3 u2(y-u)- A 4 ( --(y-u)-A 5 F 


(ku - y) 

Ae (k-1)u d(fy2) - 1/(k-1) (7.55) 

X j dy 

y'^(u-y) (u-y)y'^ dy 


which determines P, Ai, Az, A 3 , A 4 and As are constants depending on time. As 
a consequence of (7.2), (7.4) and (7.54) we write 

Bi-m 

1 1'"^^ (y 2u) I 1 1<^ (Biii^ (y (I) 1- - (U y) ■ B.r 

7 

Ku-y (k-1)u d(fy^) 

XB 4 B 5 (Ku-y )1 — dy (7.56) 

y^(u-y) (u-y)y^ dy 

where Bi, B 2 , B 3 , B 4 and B 5 are constants depending on time. From equation 
(7.2) on using (7.47) and (7.48) we get 

1 dp 1 du 2u 

- (--) = — ( + — ) (7.67) 

p dr r-u dr r 

Differentiating (7.55) with respect to r and then using (7.57) we get 
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(3-k) (r - 2u) (du/dr + 2u/r) 

p3-k ( + _ 2du/dr) + 

r-u 


(2-k) (du/dt + 2u/r) 

p2-k [ u'(r-u) + 

r-u 

□2 m (u-r) 

D2F] 

r 

d D 4 (ku-r) (k-1)u d(fr^) 

= — [ D 2 (ku-r) I dr] (7.58) 

dr r^(u-r) (u-r)r^ dr 


D3m(u-r) d 

D3F} + — (Diu^ (r-u) + 

r dr 


The equation (7.58) and (7.54) express relationship between u and r then 
(7.54) (7.55) & (7.56) express p, p and H in terms of r and therefore, given the 
required solution. 

In the absence of any radiation effects. The klimshin’s cofficient K becomes the 
usual adiabatic exponent r and then the solution referred above agrees with 
corresponding solution in this chapter. 
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